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ork  during  this  contract  has  been  directed  towards  further  im- 
proving and  evaluating  the  automated  depth  determination  procedure 
developed  during  three  previous  contracts.  This  procedure  uses 
computed  travel  times  to  combine  depth  phase  information  asso- 
ciated with  P,  PP,  PPP,  and  PcP  arrivals  recorded  at  a network  of 
stations.  The  resulting  depth  estimate  is  better  than  those  ob- 
tained by  conventional  methods,  which  do  not  make  full  use  of  all 
this  depth  phase  information. ^ 


Significant  program  improvements  inTpi-ejiifiilted  include  a set  of 
modifications  which  enhance  the  detection  of  depth  phases,  and  a 
statistical  technique  for  assessing  the  reliability  of  deoth  esti- 
mates. The  depth  determination  procedure  was  evaluated  by  apply- 
ing it  to  nine  different  events.  Results  of  this  analysis  indi- 
cate that  depth  estimates  with  an  error  of  t3  km  can  be  expected 
for  events  deeper  than  10  km  if  sesimograms  from  at  least  five 
stations  are  available. 

Finally,  the  depth  determination  procedure  has  been  implemented 
at  the  Seismic  Data  Analysis  Center  (SDAC) . This  will  permit  the 
technique  to  be  tested  on  a large  base  of  events,  for  the  purpose 
of  thoroughly  evaluating  its  place  in  the  nuclear  monitoring  pro- 
gram. 
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1.0  INTRODUCTION 


One  of  the  most  physically  straightforward  discriminants  for 
underground  explosions  is  source  depth.  If  a sufficiently 
accurate  method  of  estimating  seismic  source  depth  is  avail- 
able, those  events  deeper  than  possible  for  a man-made  source 
can  immediately  be  classified  as  earthquakes,  and  further 
analysis  can  be  concentrated  on  the  remaining  shallow  events. 

The  basic  objective  of  ENSCO's  research  in  this  area  has  been 
to  develop  a depth  estimation  technique  that  is  sufficiently 
accurate  and  apnlicable  enough  for  routine  discriminant  analy- 
sis. 

At  teleseismic  distances,  the  most  reliable  information  about 
the  depth  of  an  event  is  contained  in  the  delay  time  between 
the  depth  phases  and  their  corresponding  primary  phases.  In 
practice,  when  an  accurate  depth  estimate  is  desired,  it  is 
usually  obtained  from  the  pP-P  delay  time,  determined  by  visu- 
ally identifying  the  pP  arrival.  This  technique  works  well 
when  seismograms  with  clearly  identifiable  pP  arrivals  are 
available,  but  it  is  often  not  possible  to  pick  pP  unambiguously 
on  any  of  the  records  at  hand.  An  alternative  approach  is  to 
use  correlation  techniques  to  estimate  the  pP-P  time,  but, 
when  applied  only  to  the  first  arrival  portion  of  the  seismo- 
gram, this  method  rarely  produces  better  results  than  visual 
identification.  Our  basic  approach  has  been  to  improve  on 
these  conventional  depth  determination  methods  by  making  more 
complete  use  of  the  information  contained  in  the  depth  phases. 

In  the  three  previous  contracts,  the  basic  framework  of  an 
improved  depth  determination  procedure  was  developed.  The 
first  contract  (Project  VT/4710)  demonstrated  that  better  esti- 
mates of  p-  and  s-  phase  delay  times  could  be  obtained  by  apply- 
ing cepstrum  techniques  to  the  entire  P-wave  coda,  including 


w 


possible  PP,  PPP,  and  PcP  arrivals  and  their  corresponding 
depth  phases.  During  the  second  contract  (Project  VT/5710), 
the  cepstrum  matched  filter  (CMF)  algorithm  was  developed  to 
simplify  cepstrum  interpretation  by  combining  all  the  cepstrum 
peaks  expected  from  correlations  between  p-,  s-,  and  primary 
phases  into  a single  peak.  This  contract  also  demonstrated 
the  use  of  travel  time  information  to  correct  for  changes  in 
the  primary  phase- to-depth  phase  delay  times  for  later  phases, 
enabling  cepstrums  computed  from  different  parts  of  the  seismo- 
gram to  be  constructively  stacked.  During  the  third  contract 
(Project  VT/6710),  the  use  of  travel  time  information  was 
automated,  resulting  in  a depth  determination  procedure  that 
uses  seismograms  from  a group  of  stations  to  produce  plots  of 
stacked  cepstrum  amplitude  vs.  depth.  Ideally,  these  depth 
plots  combine  all  the  available  depth  phase  information  for 
an  event  into  a single,  easily  interpretable  display. 

The  basic  objective  of  the  present  contract  has  been  to  im- 
prove on  this  depth  determination  technique  and  to  apply  it 
to  a larger  group  of  events.  Along  with  many  smaller  program 
modifications,  two  major  additions  have  been  made  to  the  depth 
determination  procedure: 

• A new  cepstrum  peak-picking  algorithm  that 
uses  calculated  s-phase  delay  times  as  well 
as  p-phase  times  (replaces  CMF) . 

• A statistical  technique  for  estimating  the 
significance  of  depth  plot  peaks. 

The  first  of  these  new  features  is  designed  to  do  a better 
job  of  locating  the  desired  cepstrum  peaks,  and  the  second 
is  intended  to  simplify  the  interpretation  of  the  final  depth 
plots.  Concurrent  with  this  program  development,  nine  new 
events  were  analyzed,  with  source  depth  estimates  being  made 
on  eight  of  them. 
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Finally,  the  source  depth  analysis  program  has  been  imple- 
mented at  the  Seismic  Data  Analysis  Center  (SDAC) . This  is 
a totally  rewritten  program,  and  has  several  advantages  over 
previous  versions,  including: 

• easier  access  to  data 

• more  flexibility  in  varying  program  parameters 
and  analysis  options 

• more  complete  graphics  CCalcomp  plot)  output 

In  this  form,  the  depth  determination  procedure  can  be  used 
for  event  analysis  on  a larger  scale  than  has  previously  been 
practical,  thus  allowing  the  performance  of  the  procedure  to 
be  more  systematically  analyzed. 

This  report  describes  the  program  development  and  event  analy- 
sis work  done  during  the  past  year.  A brief  outline  of  the 
depth  determination  procedure  will  be  given,  with  more  detailed 
accounts  of  the  major  new  program  features.  Results  from  the 
nine  events  analyzed  will  be  presented,  along  with  their  inter- 
pretations. Finally,  the  possibilities  for  further  improve- 
ment and  application  of  this  procedure  in  the  nuclear  monitor- 
ing program  will  be  discussed. 


2.0  PROGRAM  DEVELOPMENT 
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Basically,  the  program  development  done  during  this  contract 
was  directed  towards  simplifying  the  task  of  obtaining  a source 
depth  estimate.  This  has  been  accomplished  by  improving  the 
depth  determination  procedure  in  three  ways: 
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• Reducing  spurious  depth  plot  peaks 

• Providing  a more  objective  criterion  for 
estimating  the  reliability  of  a given  peak 

• Allowing  greater  flexibility  in  changing 
key  input  parameters  and  analysis  options 

The  following  sections  describe  the  particular  program  changes 
that  have  been  made  to  obtain  the  desired  improvement  in  each 
of  these  areas. 

2 . 1 The  Seismic  Source  Depth  Determination  Procedure 

Before  discussing  particular  program  modifications,  a brief 
description  of  the  depth  determination  procedure  is  desirable; 
a more  detailed  treatment  is  given  in  the  final  report  from 
the  previous  contract. 

The  basic  components  of  this  procedure  are  outlined  in  Figure 
1.  First,  each  seismogram  to  be  processed  is  divided  into  a 
series  of  time  windows  of  equal  length,  where  the  window  length 
and  the  total  amount  of  data  to  be  used  are  user-input  para- 
meters. Next,  a cepstrum  is  computed  for  each  time  window. 
Finally,  the  travel  times  for  this  station  distance  are  com- 
puted for  each  trial  depth  and  used  to  pick  the  cepstrum  ampli- 
tudes to  be  accumulated  into  the  appropriate  individual  phase 
depth  plots. 
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The  process  is  repeated  for  each  station,  resulting  in  four 
I individual  phase  depth  plots  that  contain  depth  information 

1 from  all  the  seismograms.  When  every  station  has  been  pro- 

j cessed,  the  individual  phase  plots  are  weighted  according  to 

the  number  of  contributions  and  combined  into  a single  com- 
posite depth  plot,  which  contains  all  the  depth  information 
present  in  the  whole  set  of  seismograms.  Ideally,  the 
largest  peak  on  the  composite  depth  plot  gives  the  best  esti- 
mate of  the  source  depth. 

2 . 2 Reduction  of  Spurious  Depth  Plot  Peaks 

If  all  the  assumptions  implicit  in  the  denth  determination 
procedure  were  correct,  the  final  depth  plots  would  consist 
only  of  a single  large  peak  at  the  true  source  depth,  which 
includes  information  from  two  peaks  arising  from  correlations 
between  the  p-  and  s-  arrivals  (difference  peak)  and  between 
the  primary  and  s-  arrivals  (second  peak).  Depth  plots  gen- 
erated from  real  data,  however,  contain  many  additional  peaks 
The  occurrence  of  these  spurious  peaks  has  been  reduced  by 
making  modifications  to  both  the  cepstrum  and  depth  plot  pro- 
cessing algorithms. 

2.2.1  Cepstrum  Processing 

Efforts  were  made  to  reduce  false  cepstrum  peaks  arising  from 
two  sources:  inter-seismic  phase  correlations  and  source 
spectrum  structure.  The  particular  techniques  used  can  be 
i,< ’ • . ?d : 


• Removal  of  source  spectrum  structure 

cosine  taper  of  short  lag  part  of 
cepstrum 

spectrum  whitening 

• Reduction  of  inter-phase  correlations 


overlapping  seismogram  windows 


The  cosine  taper  resulted  in  a substantial  improvement,  and 
has  been  permanently  incorporated  into  the  depth  determina- 
tion procedure.  The  other  techniques,  however,  require  addi- 
tional investigation  to  determine  whether  they  can  be  useful. 


Cosine  Taper 


Because  of  the  large  low  frequency  content  of  the  source  wave- 
let power  spectrum,  many  high  amplitude  peaks  tend  to  occur  in 
the  short  lag  part  of  the  cepstrum,  resulting  in  spurious  depth 
plot  peaks  falling  at  shallow  depths.  One  of  these,  correspond- 
ing to  the  first  side  lobe  of  the  autocorrelation  function, 
always  occurs  at  a lag  of  1.5  - 2.0  sec.,  producing  a large 
depth  plot  peak  near  7 km.  This  peak  has  been  successfully 
eliminated  by  applying  a cosine  taper  to  the  short  lag  part 
of  each  cepstrum.  This  technique  was  tested  on  four  events: 
Illinois  (Figure  A-2),  Andreanof  Islands  (Figure  A-26,  A-28) , 
Turkey  (Figure  A-5),  and  Montana  (Figure  A-10) , with  clear 
improvement  resulting  in  all  cases.  The  success  of  this 
method  for  the  Montana  event,  with  an  estimated  depth  of  13 
km,  demonstrates  that  the  cosine  taper  does  not  degrade  actual 
depth  phase  peaks,  even  for  events  with  depths  near  10  km. 


Spectrum  Whitening 
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A more  general,  and  more  effective,  way  of  eliminating  cepstrum 
peaks  due  to  the  structure  of  the  source  wavelet  spectrum  is 
to  whiten  the  source  spectrum.  If  the  source  wavelet  had  a 
flat  amplitude  spectrum,  the  cepstrum  peaks  resulting  from 
depth  phase  correlations  would  be  delta  functions;  the  spuri- 
ous peaks  arising  from  low  frequency  modulation  in  the  source 
spectrum  would  not  be  present.  Several  whitening  techniques, 
including  log  whitening  and  inverse  filtering,  were  tried  in 
an  effort  to  eliminate  these  peaks,  but  without  success. 
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Overlapping  Seismogram  Windows 

Spurious  cepstrum  peaks  can  also  occur  if  arrivals  from  dif- 
ferent primary  phases  (e.g.,  sP-PP,  PPP-PP)  fall  into  the  same 
seismogram  window.  The  first  attempt  at  reducing  these  peaks 
involved  computing  cepstrums  using  overlapping  time  windows 
instead  of  non-overlapping  windows.  Although  this  does  not 
prevent  arrivals  from  different  primary  phases  from  falling 
within  the  same  window,  it  does  make  it  more  likely  that  there 
will  always  be  a time  window  that  contains  the  complete  set 
of  primary,  p-,  and  s-  arrivals.  On  the  small  set  of  events 
analyzed,  this  technique  failed  to  improve  on  the  results  ob- 
tained with  non-overlapping  windows. 

Arrival  Time-Controlled  Seismogram  Windows 

The  second  approach  to  the  inter-phase  correlation  problem 
used  non-uniformly  spaced,  variable  length  time  windows  to 
guarantee  that  depth  phases  would  always  fall  in  the  same 
window  as  their  corresponding  primary  phases.  For  each  pri- 
mary phase,  this  method  computes  cepstrums  from  time  windows 
starting  at  the  expected  arrival  time  for  that  phase,  and  ex- 
tending far  enough  to  include  the  corresponding  s-  arrival. 
Since  the  primary- to-s  delay  time  increases  with  source  depth, 
this  method  requires  the  simultaneous  use  of  cepstrums  computed 
from  time  windows  of  different  lengths. 

The  arrival  time  controlled  method  also  failed  to  produce  the 
desired  depth  plot  improvement.  This  technique,  however, 
still  seems  like  a good  idea.  Its  lack  of  success  may  be  due 
to  inconsistent  scaling  between  cepstrums  computed  from  time 
windows  of  different  lengths;  further  investigation  in  this 
area  would  probably  be  useful. 
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2.2.2  Depth  Plot  Generation 

False  depth  plot  peaks  can  also  be  introduced  during  the  pro- 
cedure that  generates  depth  plots  from  cepstrums.  To  minimize 
this,  several  modifications  have  been  made  to  the  depth  plot 
generation  procedure;  by  far  the  most  important  of  these  is 
the  use  of  calculated  s-phase  delay  times  in  the  cepstrum 
peak-picking  algorithm.  This  has  allowed  the  old  CMF  algo- 
rithm to  be  replaced  by  a much  simpler  one,  and  has  elimin- 
ated one  source  of  spurious  peaks. 

Basically,  the  availability  of  calculated  s-delay  times  en- 
ables the  program  to  locate  difference  and  second  peaks  dir- 
ectly, instead  of  scanning  over  a range  of  T^  ^"^p- 
was  done  by  CMF.  For  each  trial  depth  and  station,  CMF  used 
a 3-dimensional  polynomial  to  calculate  the  expected  delay 
time  for  the  main  peak,  then  scanned  over  a range  of  delay 
times  for  the  difference  and  second  peaks.  This  made  it  rela- 
tively easy  to  pick  up  high  amplitude,  spurious  cepstrum 
peaks  and  miss  the  correct  peaks  entirely.  With  the  inroduc- 
tion  of  polynomial  coefficients  for  s-phases,  however,  pre- 
calculated delay  times  can  also  be  used  for  the  difference 
and  second  peaks,  thus  reducing  the  chances  of  picking  up  a 
spurious  peak. 

2.3  Calculation  of  Significance  for  Depth  Plot  Peaks 

Depth  plots  generated  from  real  data  will,  in  general,  con- 
tain many  peaks,  and  it  is  desirable  to  have  some  quantita- 
•tive  means  of  estimating  whether  any  of  these  peaks  indicate 
the  true  source  depth.  The  distinguishing  characteristic  of 
cepstrum  peaks  that  arise  from  depth  phases  is  that  their 
locations  vary  with  station  distance  and  seismic  phase  in  the 
way  expected  by  the  cepstrum  picking  algorithm.  Thus,  these 
peaks  will  always  contribute  to  the  depth  plots  at  the  same 
depth,  i.e,  the  source  depth,  for  each  cepstrum  processed. 
Cepstrum  peaks  arising  from  other  sources  (interphase  correla- 
tions, source  spectrum  structure,  etc.)  will  not,  in  general. 
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vary  between  cepstrums  in  the  way  predicted  by  the  p-  and  s- 
delay  time  calculations.  These  false  cepstrum  peaks,  in  con- 
trast with  the  depth  phase  cepstrum  peaks,  will  not  stack 
when  summed  into  the  depth  plot  arrays.  Consequently,  the 
range  of  amplitudes  in  a depth  plot  that  contains  real  depth 
information  should  be  much  larger  than  in  a depth  plot  contain- 
ing only  false,  non-depth  phase  peaks.  Based  on  this  differ- 
ence, a quantitative  technique  for  estimating  the  signifi- 
cance of  depth  plot  peaks  has  been  developed. 

Basically,  this  technique  consists  of  using  the  actual  cep- 
strums for  an  event  to  generate  a set  of  "random  depth  plots" 
containing  no  depth  information.  The  statistical  properties 
of  these  plots  are  then  determined  and  compared  with  the  real 
depth  plots.  If  the  real  depth  plots  have  peaks  A^fith  ampli- 
tudes that  are  significantly  higher  than  those  predicted  by 
the  amplitude  distribution  in  the  random  depth  plots,  then 
those  peaks  probably  contain  true  depth  information. 

The  first  step  in  the  significance  level  calculation  is  to 
generate  a set  of  depth  plots  similar  to  those  that  would  re- 
sult from  the  cepstrums  if  they  contained  no  depth  phase  in- 
formation. This  is  done  by  choosing  the  cepstrum  values  that 
contribute  to  the  depth  plot  at  each  depth  using  random  delay 
times  instead  of  the  true  delay  times  for  the  difference,  main, 
and  second  peaks  expected  for  that  depth.  Specifically,  the 
contributon  of  a cepstrum  c(t)  to  this  random  depth  plot  R(d) 
is  determined  by  generating,  for  each  depth  d,  three  random 
numbers  t , t , x , uniformly  distributed  between  0 and  the 

l 2 3 

cepstrum  window  length.  After  arranging  the  numbers  so  that 
T <T  <T  , the  value  of  the  random  depth  plot  at  denth  d is  up- 
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dated: 


R(d)  » R(d)  + { 


c(t  )+c(t  )+c(t  ) if  0.7  max  (c(t  ),cCt  ))<c(t  ) 

12  3 13—2 

c(t  ) Otherwise 
2 


This  is  the  same  algorithm  used  to  determine  the  contribution 
of  a cepstrum  c(t)  to  the  real  depth  plots,  except  that  the 
use  of  random  numbers  nrevents  any  true  depth  phase  peaks  pre- 
sent in  c(t)  from  contributing  at  the  correct  depth.  In  other 
words,  the  use  of  random  delay  times  makes  c(t)  look  like  it 
contains  no  depth  phase  information. 

This  procedure  is  repeated  for  each  cepstrum,  until  all  cep- 
strums  have  been  processed  and  the  final  values  for  the  random 
depth  plot  R(d)  are  obtained.  This  random  depth  plot  can  be 
viewed  as  a set  of  values  of  a random  variable  r with  unknown 
probability  distribution  f(r);  i.e., 

f(r)  = Probability  R(d)  = r 

Note  that  r is  not  uniformly  distributed,  even  though  the 
random  delay  times  are.  This  is  because  each  r is  a sum  of 
cepstrum  values,  and  each  cepstrum  c(t)  is  a nonlinear  func- 
tion of  the  delay  time. 

The  next  step  is  to  use  R(d)  to  estimate  f(r)  and  then  to 
estimate  the  cumulative  distribution  function  F(r)  which  will 
be  used  to  obtain  the  desired  significance  level  values.  The 
details  of  this  process  are  illustrated  in  Figure  2. 

First,  f(r)  is  estimated  by  dividing  the  amplitude  axis  on  the 
random  depth  plot  into  intervals  of  width  Ar  and  counting  the 
number  of  points  that  fall  in  each  interval.  Then  the  esti- 
mated probability  distribution  is  given  by: 
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f (nAr) 


no.  of  points  with  amplitudes  in  ((n-l)Ar,  nAr) 
total  number  of  ooints 


Next,  this  estimate  of  f(r),  shown  in  Figure  2b,  is  used  to 
find  the  cumulative  distribution  function  F(r)  by  computing 
the  sums 

n 

F(nAr)  = I f (iAr) 
i=l 

Finally,  F(r)  is  used  to  compute  the  significance  levels  as 

shown  in  Figure  2c.  F(r)  is  the  probability  of  getting  a 

random  depth  plot  point  with  a value  <r,  and  the  significance 

levels  are  given  by  the  inverse  of  F(r).  Specifically,  the 

80%  and  95%  significance  levels  are  the  depth  plot  values  r 

and  r that  satisfy 
2 

F(r  ) = 0.80 
1 

F(r  ) = 0.95 
2 

These  two  values  are  marked  on  the  final  deoth  plot  (Figure 
2d),  giving  the  analyst  a quantitative  criterion  for  estimat- 
ing the  validity  of  depth  plot  peaks.  For  example,  if  a depth 
plot  has  a peak  amplitude  above  the  95%  level  (i.e.,  the  amp- 
litude that  is  greater  than  95%  of  the  points  on  the  random 
depth  plot) , the  peak  (or  peaks)  containing  these  points  are 
likely  to  contain  real  depth  information.  On  the  other  hand, 
a peak  that  falls  entirely  below  the  95%  level  could  be  a 
false  one,  even  if  it  is  the  dominant  peak  on  the  depth  plot, 
since  5%  of  the  points  computed  using  random  delay  times  on 
the  same  cepstrums  had  a greater  amplitude. 
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2 . 4 Implementation  at  SDAC 


After  most  of  the  program  development  work  had  been  completed, 
the  source  depth  determination  program  was  rewritten  for  im- 
plementation at  SDAC.  This  version  includes  all  the  new  fea- 
tures developed  during  this  contract,  and,  in  addition  to 
allowing  easier  access  to  data,  is  characterized  by  increased 
flexibility  in  changing  key  input  parameters  and  in  obtaining 
different  kinds  of  output.  The  basic  purpose  of  this  new  ver- 
sion is  to  provide  an  efficient  tool  for  performing  source 
depth  analysis  on  a large  number  of  events. 

The  general  floxv  of  the  SDAC  depth  determination  program  is  shown 
in  Figure  3.  The  program  consists  of  two  parts,  enclosed  by 
dashed  lines  in  the  figure,  which  can  be  executed  either  separ- 
ately or  together.  The  first  part  computes  a set  of  cepstrums 
from  the  input  seismogram  data  tape,  and  writes  these  cepstrums 
on  an  output  tape.  Control  parameters,  such  as  the  total  data 
length,  cepstrum  computation  window  length,  and  stations  used, 
are  input  on  data  cards.  Optional  Calcomp  plots  of  the  seis- 
mograms may  also  be  generated.  The  second  part  of  the  program 
reads  the  cepstrum  file  written  by  the  first  part  and  gener- 
ates a set  of  depth  plots.  The  control  parameters  used  in 
producing  the  depth  plots  may  either  be  the  ones  specified  in 
the  original  run  of  the  first  step,  or  they  may  be  new  ones 
input  to  the  second  step  on  a set  of  override  control  cards. 

This  override  feature  makes  it  possible  to  obtain  sets  of  depth 
plots  for  different  groups  of  stations,  data  lengths,  window 
lengths,  and  other  control  parameters,  without  recomputing  the 
cepstrums.  Calcomp  plots  of  individual  cepstrums  are  also 
optionally  available  from  the  depth  plot  generation  step.  The 
various  plots  output  by  both  steps  are  labelled  with  all  iden- 
tifying information  and  pertinent  control  parameters,  and  com- 
prise a complete  report  of  the  depth  analysis  process  for  each 
event . 
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3.0  ANALYSIS  OF  EVENTS  5 

i 

i 

] 

i 

Nine  events  have  been  analyzed  during  this  contract,  with 
source  depth  estimates  being  made  on  all  but  one  of  them. 

Results  of  these  analyses  are  summarized  in  Table  1;  seis- 
mograms and  key  depth  plots  for  each  event  can  be  found  in 
the  Appendix.  A reliability  grade  has  been  listed  in  Table 
1 for  each  event  interpreted.  This  grade  is  a qualitative 
estimate  of  the  reliability  of  the  corresponding  depth  esti- 
mate, and  is  based  on  factors  such  as: 


li 


I 


• appearance  of  depth  plots 

• number  of  stations  processed 

• seismogram  signal-to-noise  ratios 

• significance  level  of  peaks 

• consistency  of  peaks  among  different 
stations  and  phases 

• distribution  of  station  distances 

In  this  section,  a detailed  discussion  of  the  interpretation 
of  the  Andreanof  Islands  event  will  be  given.  This  event  was 
chosen  because  of  its  complex  interpretation,  and  because  the 
large  number  of  stations  available  allowed  a more  complete 
analysis  to  be  performed.  In  addition,  the  interpretation  of 
each  of  the  remaining  events  will  be  briefly  described. 

3.1  Andreanof  Islands  Event 

Data  for  the  Andreanof  Islands  event  consists  of  seismograms 
from  14  stations  with  A's  ranging  from  34°  to  73°;  these  seis- 
mograms are  shown  in  Figures  A-21  and  A-22.  Because  of  the 
large  number  of  stations  available,  the  data  could  be  divided 
into  smaller  groups  of  stations  and  each  group  processed 
separately.  The  resulting  depth  plots  can  then  be  used  to 
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Table  1 


Event  Analysis  Results 


NO.  OF  RELIA- 


EVENT 

STATIONS 

DEPTH 

BILITY 

FIGURES 

Illinois  (11/9/68) 

6 

26  km 

G 

A-1,  A-2 

Turkey  (3/27/75) 

4 

17  km 

G 

A-3  to  A-5 

Kuril  (3/23/75) 

3 

39  km 

F 

A-6,  A-7 

Montana  (6/30/75) 

5 

13  km 

F 

A-8  to  A-10 

Mojave  Desert  (6/1/75) 

3 

14  km 

P 

A-11,  A-12 

Idaho  (3/28/75) 

2 

20  km 

Q 

A-13,  A-14 

Iran  (3/7/75) 

2 

19  km 

Q 

A-15  to  A-17 

Kasmir  (4/28/75) 

2 

— 

-- 

A-18  to  A-20 

Andreanof  Is,  (11/22/65) 

14 

25  km 

F 

A-21  to  A-38 

1 

RELIABILITY: 

G “ Good 
F » Fair 
P « Poor 
V « Very  Poor 
Q - Questionable 


check  the  behavior  of  individual  peaks  for  different  station 
groups,  thus  providing  the  interpreter  with  more  information. 
This  technique  has  proven  to  be  quite  useful,  since  the  depth 
plots  for  this  event  contain  three  distinct  peaks  whose  valid- 
ity would  be  difficult  to  assess  from  only  one  large  group  of 
stations , 

Seven  different  station  groups  were  used  in  analyzing  the 
Andreanof  Islands  event;  the  characteristics  of  these  data 
sets  are  listed  in  Table  2.  The  resulting  depth  plots  are 
shown  in  Figures  A- 23  through  A-38.  Each  data  set  was  run 
using  the  first  102.4  seconds  of  the  seismograms,  and  depth 
plots  were  generated  using  window  lengths  of  25.6  and  51.2 
seconds . 

Four  distinct  peaks  reappear  consistently  throughout  all  the 
depth  plots,  at  depths  of  7,  15,  25,  and  42  km.  The  7 km 
peak  is  an  autocorrelation  side  lobe  that  can  be  eliminated 
by  the  cosine  taper  algorithm;  Figures  A-26  and  .4-28  (depth 
plots  from  Data  Set  2 run  with  the  cosine  taper)  show  this 
clearly.  The  other  three  peaks,  however,  appear  to  represent 
real  seismogram  features. 

Arrivals  appear  on  many  of  the  seismograms  at  delay  times 
near  7 seconds  and  10  seconds  after  P,  suggesting  that  the 
25  km  peak  should  be  interpreted  as  the  main  one,  with  the 
42  km  peak  as  the  second  peak  and  the  15  km  peak  as  the  dif- 
ference peak.  Estimates  of  the  "average"  delay  times  corres- 
ponding to  these  peaks  show  that  the  15  km  peak,  with  a delay 
time  of  about  4.5  seconds,  does  correspond  to  the  difference 
peak  between  the  25  km  (7.1  seconds)  and  42  km  (11.7  seconds) 
peaks.  However,  the  42  km  peak  is  more  than  1 second  later 
than  the  expected  sP  time  for  a 25  km  source  depth.  This 
anomalously  late  arrival  explains  why  the  cepstrum  picking 
algorithm  failed  to  sum  the  difference  and  second  peak  ampli- 
tudes into  the  main  peak. 
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Table  2.  Andreanof  Islands  Event 
Data  Sets 


DATA  SET  NO.  DHSCRIPTIOX 


5 stations,  distributed  over  entire 
available  range  of  A's 

10  stations,  distributed  over  entire 
available  range  of  i's 

4 stations,  A's  within  5°  of  A=39° 
(close  range  A's) 

6 stations,  A's  within  5°  of  A=48° 
(middle  range  A's) 

4 stations,  A's  within  5°  of  A=61° 
(long  range  A's) 

8 stations,  arrival  visible  7 sec 
after  P . 

5 stations,  arrival  visible  10  sec 
after  P 


The  long  s delay  time  is  most  straightforwardly  explained  as 
a sP  arrival  from  a source  region  with  abnormally  slow  S 
velocities.  An  alternative  interpretation  is  that  there  are 
two  good  reflectors,  such  as  water  bottom  and  water  surface, 
above  the  source,  with  a two-way  travel  time  of  about  4.5 
seconds  between  the  two  boundaries.  Then,  the  25  km  peak 
would  be  the  result  of  an  echo  off  the  deeper  reflector,  and 
the  42  km  peak  would  come  from  an  echo  off  the  shallower  re- 
flector . 

Thus,  the  interpretation  of  this  event  can  be  summarized: 

• The  15  km  peak  is  a difference  peak  arising 
from  the  25  km  and  42  km  peaks. 

• Two  interpretations  are  possible  for  the  25 
km  and  42  km  peaks: 

1.  The  25  km  peak  represents  a pP  arrival 
(and  thus  the  true  source  depth) , and 
the  42  km  peak  represents  a slow  sP 
arrival . 

2.  There  are  two  reflectors  above  the 
source,  with  the  25  km  peak  coming 
from  an  echo  off  the  deeper  one,  and 
the  42  km  peak  from  an  echo  off  the 
shallower  one. 

This  event  is  a good  one  for  illustrating  how  many  factors 
can  be  involved  in  obtaining  a depth  estimate.  Seismogram 
plots,  delay  time  information,  and  depth  plots  from  several 
different  station  groups  have  all  been  used  in  the  interpre- 
tation of  this  event.  All  this  information  should  be  easily 
available  to  an  analyst  doing  routine  depth  determinations, 
and  each  of  these  features  is  either  included  in  or  planned 
for  the  depth  determination  program  implemented  at  SDAC . 


3 . 2 Other  Events 
Illinois  Event 


The  two  events  that  were  easiest  to  interpret  and  yielded 
the  most  reliable  depths  were  Illinois  and  Turkey.  Both 
events  are  characterized  by  depth  plots  with  a single  domin- 
ant peak,  well  above  the  95%  level,  and  have  seismograms  with 
good  signal -to-noise  ratios  from  an  adequate  number  of  sta- 
tions. Several  of  the  seismograms  from  these  two  events 
show  clear  pP  arrivals  with  delay  times  that  agree  with  the 
estimated  source  depths. 

The  Illinois  event  (Figures  A-1  and  A-2)  is  the  same  one 
analyzed  during  the  previous  contract.  The  six  stations 
processed  have  source-to-receiver  distances  covering  a range 
sufficient  to  produce  significant  moveout  in  depth  phase 
delay  times.  P,  PP,  and  PPP  arrivals  are  present  in  at  least 
some  of  the  seismograms,  and  the  resulting  individual  phase 
depth  plots  all  have  their  highest  amplitude  peaks  near  26  km. 
The  depth  plot  shown  in  Figure  A-2  is  the  composite  plot  ob- 
tained using  the  cosine  taper  program  version,  and  is  the 
most  current  result  on  this  event.  The  high  amplitudes  that 
occur  between  10-15  km  are  probably  due  to  a combination  of 
difference  peaks  and  source  spectrum  structure. 

Turkey  Event 

The  Turkey  event  (Figures  A-3  through  .A-5)  has  fewer  stations 
and  a smaller  range  of  station  distances  than  Illinois,  but 
the  resulting  depth  plots  lack  the  spurious  shallow  peaks 
present  in  Illinois.  Comparison  of  depth  plots  without  and 
with  the  cosine  taper  shows  how  this  technique  reduces  the 
peaks  at  depths  less  than  10  km. 
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Kuril  Event 


The  Andreanof  Islands,  Kuril,  and  Montana  events  fall  into 
the  next  reliability  grade  below  Illinois  and  Turkey.  The 
Kuril  event  (Figures  A-6  and  A-7)  has  a well-defined  depth 
plot  peak,  and  the  three  seismograms  that  were  processed  are 
all  of  good  quality;  the  reliability  for  this  event  was  down- 
graded mainly  because  so  few  stations  were  available.  Signi- 
ficance levels  are  not  present  on  the  depth  plot  for  this 
event,  since  it  was  run  before  the  significance  level  algo- 
rithm was  implemented.  pP  and  sP  arrivals  are  visible  on 
the  seismogram  from  WH2YK,  at  delay  times  that  correspond  to 
the  estimated  39  km  depth. 

Montana  Event 

The  Montana  event  (Figures  A-8  through  A-10)  has  a main  depth 
plot  peak  that  is  not  as  clear-cut  as  those  on  the  previously 
discussed  events,  and  its  significance  level  is  somewhat 
lower.  However,  the  fact  that  five  stations  were  used  in 
the  processing  lends  more  credibility  to  this  peak.  Although 
some  of  the  seismograms  are  of  poor  quality,  the  use  of  all 
five  stations  was  still  found  to  produce  the  most  interpret- 
able result.  This  event  has  clearly  been  helped  by  the  intro- 
duction of  the  cosine  taper,  in  spite  of  its  shallow  focal  depth. 

Mojave  Desert  Event 

The  depth  estimate  for  the  Mojave  Desert  event  (Figures  A-11 
and  A-12)  is  not  as  reliable  as  the  Kuril  and  Montana  events, 
even  though  the  composite  depth  plot  shows  a good  peak,  well 
above  the  95%  level.  This  event  was  downgraded  for  three 
reasons.  First,  the  significance  level  algorithm  is  of  ques- 
tionable validity  when  applied  to  such  a small  number  of  sta- 
tions. Second,  the  three  station  distances  are  clustered 
around  25°,  so  there  is  practically  no  station  moveout  in  the 


depth  phase  delay  times.  Third,  only  two  of  the  three  seis- 
mograms have  a good  s ignal - to-noise  ratio.  Thus,  although 
the  estimated  depth  of  14  km  is  a valid  one,  it  is  not  as 
reliable  as  the  depth  plot  suggests. 

Idaho  Event 

The  Idaho  and  Iran  events  were  analyzed  using  only  two  sta- 
tions each,  and  the  resulting  depth  plots  are  cluttered  and 
ambiguous.  Depth  estimates  have  been  made  for  both  of  these 
events,  but  their  reliability  is  questionable.  The  Idaho 
event  (Figures  A-13  and  A-14)  has  no  single,  dominant  peak 
on  the  composite  depth  plot,  although  the  correct  depth 
probably  lies  somewhere  in  the  high  amplitude  region  between 
16-27  km.  Three  stations  were  available  for  this  event,  but 
only  two  were  used  in  the  final  depth  plot;  including  the 
WH2YK  seismogram  was  found  to  deteriorate  the  depth  plots. 

Iran  Event 

For  the  Iran  event  (Figures  A-5  to  A-17),  conflicting  depth 
plots  were  obtained  when  two  different  window  lengths  were 
used.  Interpretation  is  hampered  by  the  large  (>90°)  distances 
for  both  stations,  resulting  in  only  P arrivals  being  present, 
and  by  the  fact  that  the  significance  level  algorithm  is  mean- 
ingless for  such  a small  amount  of  data.  The  19  km  peak  on 
the  25.6  second  depth  plot  has  been  picked  as  the  depth  esti- 
mate, since  some  evidence  of  this  peak  can  also  be  seen  on  the 
12.8  second  plot. 

Kasmir  Event 

No  depth  estimate  has  been  make  for  the  Kasmir  event  (Figures 
A-18  through  A-20).  This  event  is  similar  to  the  Iran  event, 
but  the  17  km  peak  that  appears  to  be  the  best  depth  estimate 
on  the  12.8  second  depth  plot  is  not  present  at  all  on  the 
25.6  second  plot.  Both  depth  plots  contain  a lot  of  shallow 


depth  noise,  and  their  interpretabil ity  may  have  been  improved 
if  the  cosine  taper  algorithm  had  been  available. 
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4.0  SUMMARY,  CONCLUSIONS,  AND 

RECOMMENDATIONS  FOR  FUTURE  WORK 


4.1  SUMMARY 

Two  major  additions  have  been  made  to  the  source  depth  deter- 
mination procedure  during  this  contract:  the  calculation  and 
use  of  s-phase  travel  times  in  the  cepstrum  picking  algorithm, 
and  a statistical  technique  for  determining  the  significance 
of  depth  plot  estimates.  Several  techniques  for  reducing 
spurious  cepstrum  peaks  were  also  investigated,  and  one  of 
them,  found  to  be  very  useful  for  reducing  short  delay  time 
peaks,  has  been  permanently  implemented.  Finally,  a totally 
re-written  version  of  the  depth  determination  program  has 
been  implemented  at  SDAC,  allowing  this  method  to  be  used  for 
event  analysis  on  a larger  scale  than  has  previously  been 
practical . 

A total  of  nine  new  events  have  been  analyzed  with  depth  esti- 
mates being  obtained  on  eight  of  them.  Results  are  summarized 
in  Table  1 Cp-3-2).  Two  of  these  events  had  depth  plots  with 
a single  dominant  peak,  well  above  the  95%  significance  level, 
and  had  seismograms  with  good  signal-to-noise  ratios  from  four 
or  more  stations,  and  the  resulting  depth  estimates  were  classi- 
fied as  having  good  reliability.  Three  events  had  either  fewer 
stations  or  more  ambiguous  depth  plots  than  the  first  two,  and 
their  depth  estimates  were  classified  as  fair.  One  event  had 
only  three  seismograms,  with  one  of  them  being  very  noisy,  and 
had  a station  distribution  that  resulted  in  very  little  moveout 
in  the  depth  phase  delay  times,  so  its  depth  estimate  was  classi- 
fied as  poor.  Two  more  events  were  classified  as  questionable 
because  of  very  ambiguous  depth  plots  and  small  number  of  sta- 
tions. Finally,  no  depth  estimate  could  be  made  on  one  event. 


4.2  CONCLUSIONS 


The  results  of  the  two  principal  program  changes  made  during 
this  work  were  investigated  while  processing  the  nine  events. 
The  use  of  calculated  s-phase  delay  times  was  found  to  gener- 
ally improve  the  quality  of  the  depth  plots.  This  improvement, 
however,  is  not  large,  and  is  probably  more  due  to  elimination 
of  the  noise  introduced  by  the  time  scan  in  the  old  CMF  algo- 
rithm than  to  better  detection  of  s-phases.  The  availability 
of  s-phase  delay  times  greatly  simplifies  the  depth  plot  gener- 
ation process,  and  does  do  a better  job  of  finding  s-phases 
when  they  are  present,  so  it  is  a useful  addition  to  the  pro- 
gram. The  significance  level  algorithm  was  found  to  usually 
do  a good  job  of  estimating  the  significance  of  depth  plot 
peaks,  but  it  did  not  appear  to  work  well  when  only  small 
amounts  of  data  are  available.  This  is  to  be  expected  in  view 
of  the  statistical  nature  of  this  technique. 

The  results  from  event  analysis  also  provided  new  information 
on  the  accuracy  of  depth  estimates.  Basically,  these  events 
showed  that  the  errors  involved  in  this  deoth  determination 
procedure  are  of  two  kinds:  peak  detection  errors  and  peak 
precision  errors.  Peak  detection  error  is  concerned  ivith 
whether  the  depth  determination  procedure  has  successfully 
detected  depth  phases  and  displayed  the  depth  information  as 
an  interpretable  depth  plot  peak;  i.e.,  whether  the  analyst 
is  able  to  choose  the  correct  depth  plot  peak  for  his  depth 
estimate.  The  significance  level  algorithm  is  designed  to 
reduce  this  kind  of  error  by  providing  a quantitative  criter- 
ion for  deciding  whether  a depth  plot  peak  has  come  from  depth 
phases  or  noise.  Ideally,  peak  detection  error  should  be 
quantified  as  the  probability  of  finding  the  correct  depth 
plot  peak,  but  this  requires  a large  data  base  of  events  with 
independently  known  depths.  Thus,  a qualitative  estimate  of 
the  peak  detection  probability,  such  as  the  reliability  grades 
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given  in  Table  1,  must  be  used.  Based  on  the  events  listed 
in  Table  1,  the  probability  of  successful  peak  detection  is 
"good”  for  events  with  at  least  five  good  quality  seismo- 
grams available. 

The  second  kind  of  error,  peak  precision  error,  is  the  numeri- 
cal error  in  the  estimated  depth,  assuming  the  correct  denth 
plot  peak  has  been  picked.  This  error  depends  on  factors 
such  as  depth  phase  delay  time  errors,  seismogram  spectra, 
and  properties  of  processing  procedures  such  as  cepstrum  fil- 
tering and  power  spectrum  tapering.  For  the  set  of  events 
analyzed,  peak  precision  error  has  been  estimated  at  ±3  km. 

Combining  these  two  kinds  of  errors,  the  events  processed  to 
date  indicate  that,  for  events  deeper  than  10  km,  the  depth 
determination  procedure  will  have  a "good"  probability  of 
estimating  the  source  depth  to  ±3  km  if  at  least  five  stations 
are  available. 

4.3  REC0^P1E■\DATI0NS 

Some  additional  testing  and  program  development  is  necessary 
before  this  technique  will  be  ready  for  routine  use  in  the 
nuclear  monitoring  program.  First,  the  depth  determination 
procedure  should  be  applied  to  a large  base  of  events.  This 
will  more  firmly  establish  the  applicability  of  the  method  in 
a wide  range  of  situations,  and  will  permit  its  performance 
to  be  investigated  as  a function  of  source  parameters,  station 
distribution,  and  sesimogram  characteristics.  This  large 
scale  testing  is  the  function  of  the  program  version  that  has 
been  implemented  at  SDAC.  Next,  the  possibilities  for  further 
program  development  discovered  during  this  contract  need  to 
be  investigated.  These  include  additional  work  on  the  signi- 
ficance level  algorithm  and  the  use  of  phase  information  in 
the  depth  phase  detection  procedure. 
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pr 


Finally,  the  analysis  of  the  Andreanof  Islands  event  has  shown 
that  the  depth  determination  procedure  should  eventually  be 
put  in  a fully  interactive  form.  Routine  processing  of  seis- 
mic events  by  an  analyst  requires  easy  access  and  display  of 
information  such  as  sesimograms,  travel  times,  and  results 
from  individual  stations,  as  well  as  the  capability  of  obtain- 
ing rapid  output  using  new  program  parameters.  All  these  fea- 
tures will  be  available  if  the  depth  determination  procedure 
is  implemented  on  an  interactive  graphics  system. 


Illinois  Event  Seismograms 
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Turkey  Event  (3/27/75),  Composite  Depth  Plot 
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Turkey  Event  Composite  Depth  Plot,  Cosine  Taper 
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Kuril  Event  (3/23/75)  Seismograms 


Montana  Event  (6/30/75)  Seismograms 


Iran  Event  (3/7/75)  Seismograms 


Iran  Event  (3/7/75),  P Depth  Plot 


Iran  Event  (3/7/75),  P Depth  Plot 


Cepstrum  window  length  = 25. 
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Figure  A20 

Kasmir  Event  (4/28/75),  P Depth  Plot 
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Figure  A21 

Andreanof  Islands  Event  Seismograms 
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Composite  Depth  Plot,  Data  Set  2,  Cosine  Taper 
Andreanof  Islands  Event 
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Composite  Depth  Plot,  Data  Set 
Andreanof  Islands  Event 


CEPSTRUM  WINDOW  LENGTH  = 25.6  SEC 
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Andreanof  Islands  Event 


CEPSTRUM  WINDOW  LENfiTH  = 51.2  SEC 
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CEPSTRUM  WINDOW  LENGTH 


Composite  Depth  Plot,  Data  Set 
Andreanof  Islands  Event 


Composite  Depth  Plot,  Dita  Set 
Andreanof  Islands  Event 
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